UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 

LIBRARY ROUTINE K 11 - 192 



TITLE Quartimax Orthogonal Rotation of Factors (DOI Only) 

DURATION .36 mp + ct = time in seconds 

The first term is the time required for input and out- 
put where m is the number of factors, and p the number of variables in 
the matrix "being rotated. The second term is the time required for 
computation. Only very rough estimates can be made of the number c 
of cycles required for convergence. As the number of factors increases, 
with a fixed number of variables, £ increases rapidly. Centroid factors 
will^usually require a larger £ than principal axes factors. Approximate 
estimates of c may be taken from the following table: 



NUMBER of FACTORS 



k 6 8 12 16 



k 7 12 21 56 

The time t required for each cycle depends primarily on the number of 
factors and secondarily on the number of variables. The table below 
can be used for a rough estimation of t, 



NUMBER of FACTORS 


3 


5 


9 


13 


IT 


21 


t (in seconds) 


1 


IT 


50 


120 


200 


320 



DESCRIPTION The aim of rotation may be taken as that of decreasing 
the complexity of the factorial description of the variables . The 
quartimax method of rotation provides for finding the orthogonal 
transformation which maximizes the variance of the squared factor 
loadings (thereby achieving a high degree of inequality among the 
factor loadings ) . This can be shown to be equivalent to maximizing 
the fourth powers of the loadings. In applying the method, two factors 
are operated upon at a time. The maximizing angle of transformation 
for two factors i and j is given by the equation: 



- 2 - 

k \ f ki f kJ (f ki " f kj* ^ 
= l/4 arc tan 



^( f k i 2 - f >i 2 ) 2 -^ k i 2f i 2 



j 



"kj ki kj 

where f . is the matrix of factor loadings. This program obtains 

the quart imax solution "by successive transformations of pairs of 
factors until the sum of fourth powers no longer increases . For a 
more detailed account of the method ., see! J. 0. Neuhaus and 
C. F. Wrigley* "The quartimax methods An approach to orthogonal 
simple structure," Brit * J. statist . Psychol * , Nov. 195^* 
OUTPUT The output consists of a two-hole delay character 

punched at the end of each cycle of rotations and the set of rotated 
factors punched "by columns with an N at the end of each factor. 

TAPES USED 1. Program Tape. This is the same for all problems 

and may be copied from the library tape. 

2 . Parameter Tape . This is punched as follows : 

00 6K 

00 (p)F 00 (p)F 
00 27K 

00 (m)P 00 (m)F 
22 278N 
where (p) is the number of variables and (m) 
is the number of factors . 

3. Matrix Tape. Multiply each element of the matrix 
of factor loadings by one-tenth and then punch the 
scaled matrix factor by factor in the form required 
for Illinois Library Routine N 3* The character N 

is punched after the last loading on each factor. 

INSTRUCTIONS 1. Program tape 

to ILLIAC „ -. . - T, n 

^.—^.m^, 2. Parameter tape Bl. up 

OPERATOR * 

3. Data tape Bl. down 

If continuous punching has not begun when estimated 
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time has elapsed, raise the black switch to OBEY. Wait for Illiac 
to stop on the order 2k 090$ Then start by raising white switch. 
Continuous punching will begin immediately and last less than k 
minutes . 

Program ends on OF from ll6. 

CAPACITY (m) ' (p) < 7^5 

COMCENTS At the end of each complete cycle of rotations the 

increase in the sum of fourth powers, Z. E. f . ., is inspected to 

see if another cycle of rotations should be made. If the increase 

-lt-5 it- 
has been less than T x 2 x 10 then no further iterations will 

be made, where T is an integer less than 1,000. The tolerance, T, 

is taken to be 10 unless specified. If the tolerance desired is not 

ten then insert the following immediately before the "22 278N" on 

the parameter tape: 

00 291c 

00 F 00 (T) F 
If it is necessary to complete the rotation during a second time 
period on the computer, use the following procedure. The black 
switch is raised to the H obey" position about five minutes before 
the first time period is to end; the computer will complete a cycle 
of rotations and then stop. Start with the white switch as soon 
as the computer stops and the rotated factors will be printed. Proceed 
in the second session as in the first, but using the rotated factor 
tape as the matrix tape. 

Kll contains an old obsolete version of T2 which does not agree with 
the serial copy of T2 S This old version is correct^ however 

RT; IO/13/6O 



c 
DAT E November 17, 1955 

PROGRAMMED BY J. 0. Neuhaus 
CHECKED BY 



J0K. 

RJT* 



inge 



APPROVED BY 




This program was prepared for 
the library by the Psyeholojpr 
Department . 



USE OF STORAGE - QUARTIMAX 







LOCATION 


CONTENTS 


5 


00 IF 00 IF 


6 


00 pF 00 pF J variables (on parameter tape) 


7 


kO ( 2<79 ) F 00 F; cos -£- 
v +mp o 


8 


jt 
sm — 7T- 


9 


if 

" sln T 


10 


V> k 

previous L f . . « 



12 

13 
Ik 

15 
16 

17 
18 

19 

20, 

21 

22 
23 

2k 

25 

26 



i,J 1J 



-59 vi 4 

11 279 x 2 ^ during input*,, last L f. .. 



during calculation i j 



counter 
i 


(accumulati: 


ng var: 


Lables ) 


binary j 


switch 


(after forming 


sum) 


*\ 4 












^4 












k ki 


kj 










^\i 


f 5 










50 mp 


00 F 




(^ 






temporary stoi 


-e for 


f 
ki 


f *j 


Denominator 










binary £ 


switch 










l/2 cos 


jt 










l/2 sin 


it 










it 













merator 



Use of Storage - Quartimax (continued) Page 2 

LOCATION CONTENTS 

27 ni o factors (on parameter tape) 

28 50 2T9F 75 279F 

29 tolerance T 

31 DNI Library Routine N - 3 

52 SCP Library Routine P - 6 

65 sin - cos Library Routine T - 1 

95 arctan-Library Routine T - 2 

124 body - !<fe in Routine 



PLOW CBftRT - QUAHT'JMX 



f oi*b dependent 
parameters 15^-l6l 



Input factor 
matrix 0~k \ 



—HI 



and print 4-13 



test for cvgce 



t [(old sum + i+foqleranee 
Form E, , f x 



tolerance - neWiFL— 



i,J iJ laum) 14-17 I— Hrii^^l 



^•*m 'W M>M W .»iom' i i i |iT i , i, ii j jj i fju— i — i 



increase was 
less tins 



J 



print ."''"^ 



rotated 
factors 
143-153 



at. i Qif'j?* 

r ^154 



ft -*9 22, 23 



reset addresses for 
a new cycle 134™l42 



* 



reset 120*125 " 
accuisulator and 
rotation addresses 

for $ factor 





nil ..» jiWM 

... *«Uiii«mrri r 'in. 



■■m>L . ,«lp 



J 



*miim**mmmrmMmmmmmm 



set switch 

£2 



jtest switch 107,108 



reset addresses for : 3 ~» j + 

A .th . . 328-133 I 118 

i factor ,^» ,^„ ,__..— .v 

124-125 



clear accumulators 
103-107 




accumulate products 

and store at 16-19 

23 - 39 



"Wwiwiimwihimii 9\mkmn 



rotate factors 1 
i an d J 87-104 



form and store 
Numerator as H (21 ) 
or Denominator as N (22) 
____J*0~43 



g-r* *— "». .-i imm i cmB^fci ' » ' .^ . 



transfer 



; of I 45-48 



set 



cos, sin, -sin ^| switch! have Just 



calculated! 
! numerator 



mmmmmmmmmmim* » <>m 



test switch 
45 



clear switch, 
at 23 65 

* ' — 



plant cos, 
sin, -sin at 
7, 8, 9 62-65 



„-__ find s: 
+ 60-62 



1/4 arctan 
7^, 75 




■JS"-/*^ ' " """■ ' """""I test""""^'" Q-'j fj have Just calculated] 

2 \V\ r^-P-^ - — lt^ aen^lnat 9 r__J 

'■u^l^: 1 . I T — -- jtest lD |fr-|N| 49,50 



' test N = D 
51 



H (o) -* o 
53 



I find l/4 arctaii "~ 

and store as H (o) 55-57 



LOCATION 


ORDER 


Code XI 


00 


UK 


00 


F 


00 


279F 


00 


3K 


00 


F 


00 


k¥ 


00 29K 


00 


F 


00 


10F 


00 


5K 


00 


IF 


00 


IF 


00 


7K 


ko 


SS 


00 


F 


00 


20K 


50 


SS 


00 


F 


00 


2kYL 


OOF00461939766200J 


00F001913 1 +17l6l50J 


OOF00592699681699J 


00 


28K 


50 


SS 


75 


SS 


00 


31K 


Library 


Routine N 3 


00 


52K 


Library 


Routine P 6 



NOTES 



Decimal order input 



PAGE 1 



Kil 



-Print parameter 



■tolerance 



1/2 cas — 
8 

1/2 sin ■ 

8 

8 



Decimal number sequence input 



Single column print 



LOCATION ORDER 


nomas mge 2 


00 65K 
Library Routine T 1 


Sine, cosine 


00 95K 




Library Routine T 2 


Arctan 




00 124K 







ko ss 




.■ 


50 L 




1 


26 31F 


. 




L5 L 




2 


L4 6f 
46: L 


Input factor matrix 


1 

3 


10 ?F 






36 % 




4 


26 L 






4i igf 




5 


50 ss 




■ ! 


75 ss 




6 


00 2F 


; 




40 F 




T 


50 F 

75 F 




8 


00 2F 






L4 10F 




9 


40 10F 
Up 5L 


4 
Form sum of fourth powers (E. ,f . . ) 

■■ 
test for convergence, store 


10 


L4 5F 
40 5L 


11 


LO 20F 
32 12L 


as N(10) 


12 


26 5L 






92 515F 




13 


L5 11F 






L4 29F 




^ 


LO 10F 






32 16L 


'■ ! 



Kll 



LOCATION 

16 

17 
18 

19 
20 
21 
22 

23 
2k 

25 
26 
27 
28 
29 
30 



ORDER 



kw 



L5 10F 
kO 11F 
2k 20L 
92 59F 
92 59F 
92 135F 
26 143L 
00 F 
00 F 
00 F 
kl 12F 
L5 5F 
40 13F 
kO l^F 
k9 15F 
26 23L 

50 ss 
75 SS 
00 IF 
kO F 
iA i6f 
kO 16F 
50 SS 
75 SS 
00 IF 
kO IF 
Lk 17F 
40 17F 
50 SS 
75 SS 
40 2F 
50 F 



NOTES 
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Kll 



Transfer to print routine if started 
■with white switch 



k -*0 



i, j -»1 



increase 2 E, f.. = N (l6) 



2 
increase 2 L^ f, . = N"(l7) 



MMHWaMMMMMM 



LOCATION 



31 
32 
33 
34 
35 
36 

-^ i 

38 

39 
40 
41 
42 

43 
44 
45 
46 
47 



ORDER 



7J 2F 
L4 18F 
40 18F 
50 IF 
7J 2F 
L4 19F 
40 19F 
L5 23L 
L4 5F 
40 23L 
46 29L 
L5 26L 
l4 5F 
4o 26l 

42 29L 
L5 12F 
L4 5F 
40 12F 
LO 6f 
36 42L 
26 23L 
00 F 
L5 18F 
LO 19F 

40 21F 
LI 15F 
36 4% 
49 23F 

41 15F 
L5 24F 
40 7F 
L5 25F 
40 8F 
LI 25F 



NOTES 
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KL1 



increase 4 T^ f^ ± f = N(l8) 



increase 4 Z. f,. f~ . = N(l9) 



Accumulate 

Sums 



advance addresses 



test for completion of sums 



store Numerator as N(l9) 



test switch 



set switches 



store cos, sin, - sin of —rr- 
1 at 1, 8, 9 



LOCATION 


ORDER 


NOTES PAGE 5 


kQ 


kO 9F 


-j 




26 87L 






±9 


L7 21F 


test ( | Denon . | - | Numer . | ) 




L2 22F 






50 


56 51L 
26 55L 






51 


^0 F 
LI F 


test (N = D) 




52 


52 82L 
26 77L 






55 


hi F 








50 F 






5^ 


LI 21f r 
66 22F 


(-N) -7- D 




55 


S5 F 
50 55L 






56 


26 95F 
10 2F 


find 1/4 arc tan 




57 


kO F 
L5 21F 


test N 




58 


56 72L 
L5 22F 


test D 




59 


56 68L 


test D - N 


Find the angle 




LO 21F 




of rotation 6 


60 


52 66l 

L5 F 






6l 


LO 26F 
50 6lL 


find sin [ ] 




62 


26 65F 
L5 2F 






65 


kO 7F 








L5 IF 







KL1 



LOCATION 


ORDER 




NOTES PAGE 6 


6k 


^0 8F 








LI IF 


store cos, sin, - sin 




65 


ho 9F 
hi 25F 






66 


2.6 87L 
L5 F 






61 


LO 26F 
26 6lL 






68 


iA- 2 IF 

56 71L 






69 


L5 26F 








L^ 26F 




Find the angle 


70 


L4 F 
22 6lL 




of rotation 6 


71 


25 69L 
00 F 






72 


L5 22F 
52 75L 






75 


\X 21F 
52 7^L 






7^ 


22 60L 
L5 F 




■ 


75 


22 6lL 
LO 21F 






76 


56 71L 






77 


22 7hL 
hi F 
50 F 


_« 




78 


LI 22F 
66 21F 






7Q 

1 y 


S5 F 








50 79L 







KL1 



LOCATION 


ORDER 


NOTES 
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80 


26 95F 
10 2F. 






81 


kO F 
LI F 






82 


26 57L 
50 21F 






83 


7J 22F 
32 85L 






8^4- 


LI 26F 
10 IF 






85 


26 5TL 








L5 26f 


— * 


86 


10 IF 
26 57L 






87 


hi 2F 
50 SS 






88 


7J 7F 

40 F 






89 


50 ss 
7 J 8f 






90 


l4 f 

40 IF 






91 


50 SS 
7J 9F 






92 


40 F 








50 SS 


find k component 




93 


7J IF 
L^ f 


of the rotated vectors 




9h 


40 3F 




Rotate the i 




L4 5F 




and j factors 


95 


40 SS 
L5 IF 







KL1 



LOCATION 


ORDER 


NOTES PAGE 8 


96 


L4 IF 
4o SS 






97 


F5 87L 
42 87L 






98 


42 96L 
00 20F 






99 


46 91L 
F5 92L 






100 


42 92L 








00 20F 


reset addresses 




101 


46 89L 
46 95L 






102 


L5 5F 
L4 2F 






103 


40 2F 
LO 6f 






104 


36 105L 






105 


22 87L 
4l i6f 
ko 17F 


— 


r 


106 


40 18F 








ko 19F 


clear accumulators 




107 


40 12F 








LI 23F 


test switch 




108 


B.6 118L 
L5 43L 






109 


Lk 5F 
46 43L 






110 


L5 23L 
LO 6f 






ill 


kO 23L 
46 29L 







KL1 



LOCATION 


ORDER 


112 


42 QJL 




42 S&L 


113 


46 91L 




L5 26L 


111* 


LO 6F 




40 26L 


115 


42 29L 




42 92L 


116 


46 8% 




46 95L 


117 


26 23L 




00 F 


118 


L5 l4F 




L4 5F 


119 


40 l4F 




LO 27F 


120 


36 126L 




L5.23L 


121 


LO 6f 




40 23L 


122 


46 29L 




42 87L 


123 


42 96L 




46 91L 


124 


L5 ^3L 




LO 5F 


125 


46 4JL 




26 22L 


126 


L5 13F 




L4 5F 


127 


40 13F 




LO 27F 



NOTES 



reset addresses in 



accumulation and 



rotation routines 
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KL1 



test (i - m) 



reset addresses in 



accum. and rotation 



th 
routines for j factor 



i -* i + 1 



test (i - m) 



Reset addresses, 
increase counters, 
and test to find 
position in pro- 
gram 



LOCATION 
128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 
l40 
l4l 
142 
143 



ORDER 

36 13^L 
L5 13F 
L4 5F 
40 l4F 
L5 23L 
L4 6f 
40 26L 
42 29t 
42 92L 
46 89L 
46 95L 
26 124L 
L5 28F 
40 5L 
40 23L 
46 29L 
42 8.7L 
42 96L 

46 91L 
L4 6f 
40 26L 
42 29L 
42 92L 
46 89L 
46 95L 
L5 ^3L 
LO 5F 
46 43L 
22 4l 
00 F 
4l 7F 
41 8f 



NOTES 
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KL1 



reset addresses for 



i th factor 



reset addresses to 



prepare for a new cycle 



LOCATION 



ORDER 



Ikk 


L5 SS 




50 iHl 


1^5 


26 52F 




L5 3ML 


lk6 


Lk 5F 




46 l44L 


ikj 


L5 8f 




L4 5F 


ika 


4o 8f 




lo 6f 


1^9 


36 150L 




26 ]Ml 


150 


92 7TOF 




92 129F 


151 


L5 7F 




l4 5F 


152 


kO 7F 




LO 27F 


153 


36 l^kh 




22 143L 


154 


OF F 




kl 11F 


155 


50 6f 




75 27F 


156 


00 19F 




Lk 20F 


157 


$6 20F 




46 7F 


158 


L5 23L 




L4 6f 


159 


40 26L 




k2 29L 



NOTES 
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KL1 



■ Punch 
Rotated 
Factors 



form mp 



and plant 



plant address of 1st 



Form 
dependent 

parameters 



LOCATION 



160 
161 



ORDER 



k2 921* 
k6 89L 
kG 95L 
2k L 

Sum check 



2k 999N 



NOTES 
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KL1 



loading in 2nd factor 



JON 

Rjr mge 

11/17/55 



